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Q_ij Abstract 

We find an analytic solution of the Bethe Ansatz equations (BAE) for the special 
case of a finite XXZ spin chain with free boundary conditions and with a complex 
surface field which provides for U q (sl(2)) symmetry of the Hamiltonian. More precisely, 
r> | we find one nontrivial solution, corresponding to the ground state of the system with 

anisotropy parameter A = i corresponding to = —1. 

Dedicated to Rodney Baxter 
on the occasion of his 60th birhday. 

It is widely accepted that the Bethe Ansatz equations for an integrable quantum spin 
chain can be solved analytically only in the thermodynamic limit or for a small number 
of spin waves or short chains. In this letter, however, we have managed to find a special 
solution of the BAE for a spin chain of arbitrary length N with N/2 spin waves. 

It is well known (see, for example |]] and references therein) that there is a correspondence 
between the Q-state Potts Models and the Ice- Type Models with anisotropy parameter 
A = The coincidence in the spectrum of an N-site self-dual Q-state quantum Potts 
chain with free ends with a part of the spectrum of the U q (sl(2)) symmetrical 2N-site XXZ 
Hamiltonian ([]]) is to some extent a manifestation of this correspondence. 



H xxz = Wn^n+l + °n °n+l H 7 «+l H 7 « ~ <+l)}> 

n=l 4 4 
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where A = (q + q~ 1 )/2. This Hamiltonian was considered by Alcaraz et al. [|1| and its 
U q (sl(2)) symmetry was described by Pasqier and Saleur ||. The family of commuting 
transfer-matrices that commute with H xxz was constructed by Sklyanin || incorporating a 
method of Cherednik 0]. 

Baxter's T-Q equation for the case under consideration can be written as || 

t{u)Q{u) = (j)(u + rj/2)Q(u -7]) + (j>{u - rj/2)Q{u + rj) (2) 

where q = expir/, (f)(u) = sin 2m sm 2N u and t(u) = sm2uT(u). The Q(u) are eigenvalues of 
Baxter's auxilary matrix Q(u), where Q(u) commutes with the transfer matrix T(u). The 
eigenvalue Q{u) corresponding to an eigenvector with M = N/2 — S z reversed spins has the 
form 

M 

Q( u ) = II sin ( M _ U m) sin ( M + u m)- (3) 
m=l 

Equation (|]) is equivalent to the Bethe Ansatz equations 

sin(n fc + ry/2) 
_sin(w fc - r)/2) 

In a recent article || it was argued that the criteria for the above mentioned correspon- 
dence is the existence of a second trigonometric solution for Baxter's T-Q equation and it 
was shown that in the case rj = tt/4 the spectrum of H xxz contains the spectrum of the Ising 
model. In this article we limit ourselves to the case 77 = it/ 3. This case is in some sense triv- 
ial since for this value of 77, H xxz corresponds to the 1-state Potts Model. We find only one 
eigenvalue T (u) of the transfer- matrices T(u) when Baxter's equation (Q) has two indepen- 
dent trigonometric solutions. Solving for T(u) = T (u) analytically we find a trigonometric 
polynomial Qo(u) the zeros of which satisfy the Bethe Ansatz equations (f|). The number of 
spin waves is equal to M = N/2. The corresponding eigenstate is the groundstate of H xxz 
with eigenvalue E — |(1 — N), as discovered by Alcaraz et al. |1[ numerically. 

When does a second independent periodic solution exist? This question was considered 
in article M. Here we use a variation more convenient for our goal. 

Let us consider T-Q equation (H) for 77 = f;, where L > 3 is an integer. Let us fix a 
sequence of spectral parameter values v k = vq + r)k, where k are integers and write <pk = 
(j>(yk — 77/2), Qk = Q{vk) and tk = t(vk). The functions 4>(u), Q{u) and t{u) are periodic 
with period tt. So the sequences we have introduced are also periodic with period L, i.e., 

fc +L = 4>k, etc.. 

Setting u = Vk in (||) gives the linear system 

tkQk — 4>k+xQk-i + (pkQk+l- (5) 

The matrix of coefficients for this system has a tridiagonal form. Taking vq 7^ where m 
is an integer, we have <frk 7^ for all k. 

It is straightforward to calculate the determinant of the L — 2 x L — 2 minor obtained 
by deleting the two left most columns and two lower most rows. It is equal to the product 
—0^0203 • • • 0l-i? which is nonzero, hence the rank of M cannot be less than L — 2. Here 



tt sin(K fc - n m + 77) sin(w fc + u m + 77) 

■fA sin(u fc -u m -rj) sin(w fc + u m - 77) ' 



(4) 
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we are interested in the case when the rank of M is precisely L — 2 and we have two linearly 
independent solutions for equation (|J). Let us consider the three simplest cases L = 3, 4 
and 5. The parameter 77 is equal to ~, | and | respectively. 

For L = 3 the rank of M is unity and we immediately get to = — 02, t\ = — 0o and 
£2 = Returning to the functional form, we can write 

T Q {u) = t (u)/ sin 2u = — 0(u + 7r/2)/ sin 2m = cos 2Ar u. (6) 

This is the unique eigenvalue of the transfer-matrix for which the T-Q equation has two 
independent periodic solutions. It is well known (see, for example, ||) that the eigenvalues 
of H xxz are related to the eigenvalues t(u) by 

t'(ri/2) 

E = - cos ri(N + 2- tan 2 77) + sin ri , , x . (7) 

£(77/2) 

For the eigenstate corresponding to eigenvalue (Q) we obtain E = 3/2(1 — N). This is the 
groundstate energy which was discovered by Alcaraz et al. numerically. 

Below we find all solutions of Baxter's T-Q equation corresponding to T{u) = T Q (u). 
Zeros of these solutions satisfy the BAE (f|). In particular we find Qo(x) corresponding to 
physical Bethe state. 

For L — 4, deleting the second row and the forth column of M we obtain a minor with 
determinant — 0o03(£o + h)- It is zero when £ 2 = —to, i-e., t(u + f ) = —t(u). Considering 
the other minors we obtain the functional equation 

t(u + tt/8)£(m - tt/8) = (f)(u + vt/4)0(m - tt/4) - 4>{u)(f>{u + vr/2). (8) 

This functional equation was used in || to find t(u) and show that this part of the spectrum 
of H xxz coincides with the Ising model. It would be interesting to find a corresponding Q{u). 

Lastly for L = 5, minor M35 (the third row and the fifth column are deleted) has deter- 
minant 0o04 (toti + 01^3 — 0o02)- Setting this to zero we have 

t(u)t(u + tt/5) + 0(w + 7r/10)t(u + 3tt/5) - 0(n - tt/1O)0(u + 3tt/10) = 0. (9) 

It is not difficult to check that in this case all 4 x 4 minors have zero determinant and that the 
rank of M is 3. Thus we have two independent periodic solutions of Baxter's T-Q equation. 

Note that this functional relation coincides with the Baxter-Pearce relation for the hard 
hexagon model [fT0[| . We have obtained the same truncated functional relations that have 
been obtained in |§ with the same assumptions. 

We now consider the solution of Baxter's Equation for 77 = ^ and T = T . For 77 = | and 
transfer-matrix eigenvalue T (u) = cos 2N u, the T-Q equation (|2|) reduces to 

<f>(u + 3r]/2)Q(u) + <f>(u - rj/2)Q(u + rj) + <P(u + rj/2)Q{u -rj) = 0. (10) 

This equation can be rewritten as 

f(v) + f(v + 2vr/3) + f(v + 4tt/3) = 0, (11) 
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where f(v) = sinv cos 2N (v/2) Q(v/2) has period 2ir. The trigonometric polynomial f(v) is 
an odd function, so it can be written 



K 



f( v ) = J2 c k smkv, (12) 

k=l 

where K is the degree of fiy). Then equation ( |TT|) is equivalent to c^ m = 0, m e Z. 
The condition that f(v) be divisible by sinv cos 2N (v /2) is equivalent to 

(^)7(«)Ur = 0, z = 0,l,. ..,2JV. (13) 
For even i this condition is immediate, whereas for i = 2j — 1 we use (0) to obtain 

S (-l)*^- 1 = 0, j = l,2,...,iV. (14) 

fc=l,fc^3m 

Our problem is thus to find {c^} satisfying the last equation. This problem is a special 
case of a more general problem which can be formulated as follows. Given a set of different 
complex numbers X = {xi, x 2 , . . . , xj} we seek another complex set B = /3 2 , . . . , Pi} 
where ft ^ for some i, so that 

^ftP(x,)=0 (15) 

i=l 

for any polynomial P(x) of degree not more than N — 1. It is clear that for I < N the 
system B does not exist. If (3\ ^ 0, for example, the product (x — x 2 )(x — x 3 ) . . . (x — xj) 
provides a counterexample. 

Let I — N + 1. We try the polynomials 

N 

P r = II r = l,2,...,iV. (16) 

Condition dl5|) gives f3 r P r (x r ) + PjP r (xj) = and we immediately obtain 

iV+l 

ft = const (x r — Xj) -1 , (17) 

which is a solution because the system ( |i~6"D forms a basis of the linear space of iV — 1 degree 
polynomials. So for / = iV + 1 we have a unique solution (up to an arbitrary nonzero 
constant) given by (|TTD. It is easy to show that for / = N + v we obtain a //-dimensional 
linear space of solutions. 

Returning to QT4] ) we consider = 2n, n a positive integer. Fix I = N + l = 2n + l. 
The degree K becomes 3n + 1. It is convenient to use a new index k = \3k + 1\ , where 
\k\ < n. Equation (14) can be rewritten as 

Y, ^(3« + l) 2{j - l) = 0, j = 1, 2, . . . , N, (18) 
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where we use new unknowns (3 K = (— 1) k C| 3k+ i||3k + 1| instead of Ck- Using and (|12|) we 
obtain the function f{v) 



2n + | \ / 2n 2 



/(«) = £(-!)"(, n l:){ ; + K 3 jsin(3 K + l>. (19) 

We recall that the solution of Baxter's T-Q equation for T(u) = T (u) is given by 

Q ( u ) = f (2m) / (sin 2u cos 2N u) (20) 

and its zeros satisfy the BAE (|4|). 

Another way to derive the above solution is to observe that the function f(v) satisfies a 
simple second order linear differential equation. Indeed, it is easily seen that the functions 
F + (x) and F~(x), where 

F+(x) = ± (-If ( 2 ™ + j ) ( ^ ~ i and F~(x) = F + (l/x). (21) 

are the two linearly independent solutions of the differential equation 

{((0 + n) 2 - l/9)/x + (6- nf - 1/9}F+ = 0, (22) 

where 6 = Now the fact that there is a combination f(v) of F + (e 3lv ) and F~(e 3lv ) 

which vanishes to order 2N + 1 at v = 7r follows immediately from the fact that x = — 1 is 
a singular point of the differential equation fl22"|) and that the indicial equation at this point 
has roots and 2n + 1. In terms of the variable t> , equation ( p2"D becomes 

-4 + 6n tan(3w/2) -f - + (1 - 9n 2 ) f = 0. (23) 

The zeros of f{v), the density of which is important in the thermodynamic limit, are 
located on the imaginary axis in the complex f-plane. So it is convenient to make the change 
of variable v = is. It is also useful to introduce another function g(s) = f(is)/ cosh 2n (3s/2). 
The differential equation for g(s) is then 

Let g(so) = 0. For large n we have in a small vicinity of So an approximate equation g'^u^g = 
0. This equation describes a harmonic oscillator with frequency Uq = 3n/ cosh(3so/2). The 
distance between nearest zeros is approximately As = tt/uj and we obtain the following 
density function which describes the number of zeros per unit length 

p(s) = 1/As = cu/n = 3n/(n cosh(3s/2)). (25) 



^^Up to a change of variables this is just the standard hypergeomctric differential equation, and in fact 
F+{x) = const F(-2n, 2/3 - 2n,5/3, -^a; 1 / 3 "" 
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We note that equation ( f24"|) has a history as rich as the BAE. Eckart (TT| used the 
Schrodinger equation with bell-shaped potential V(r) = —Gj cosh r for phenomenological 
studies in atomic and molecular physics. Later it was used in chemistry, biophysics and 
astrophysics, just to name a few. For more recent references see, for example, fl~2[ . 
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